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Abstract 

The superconductive phase transition (or Coulomb-Higgs phase transition 
of scalar QED in 3D) is discussed in a dual formulation which focuses on the 
magnetic rather than the electric excitations of the system. Renormalization 
group analysis of the dual formulation reveals the transition to be of second 
order and of a new universality class. Whereas coherence length and specific 
heat have XY-model exponents, the magnetic penetration depth shows mean- 
field behavior. Experimental evidence for these predictions is discussed. 
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Figure 1: Effective potential up to one-loop order. 

The study of Theoretical studies of the superconductive phase transition in the 
Ginzburg-Landau theory (or Coulomb-Higgs phase transition of scalar QED in 3D) 
have a long history. In the the early seventies, Halperin, Lubensky, and Ma 
analyzed the transition in a on-loop e-expansion and noted the absence of an infrared 
stable (IR) fixed point, suggesting that the transition be first order. A similar 
conclusion had been reached earlier in 4D by Coleman and Weinberg ||. 

The Lagrange density describing scalar QED reads, in the Euclidean formulation, 

C = \{d»A v - d u A,f + \{d, - iqA,)cft\ 2 + mj|0| 2 + A|0| 4 , (1) 

where is the electromagnetic gauge field, and (ft the self-interacting order or Higgs 
field with charge q > 0, mass and coupling constant A. As in superconductivity, 
the mass term is assumed to depend on some external parameter such as temperature 
T, changing sign at some critical value T c . Thus we set = £,q 2 t = ^ 2 (T/T C — 1), 
where £ is the length scale of the amplitude fluctuations of the Higgs field. 

At the mean- field level, the theory undergoes a transition at T = T c . The high- 
temperature Coulomb phase is characterized by a zero average of the Higgs field, 
while in the low-temperature Higgs phase this average is non-zero. For definiteness 
we consider the London limit, where the Higgs field is assumed to have a constant 
amplitude <p(x) = -j=(ft exp[iQ(x)]. The integration over the phase field 6 is easily 
carried out and results in a mass term ^m^A 2 for the gauge field, with m 2 A = q 2 (ftl, 
if we work in the Lorentz gauge d^A^ = 0. If we integrate out also the gauge field 
in the partition function, we obtain a one-loop effective potential 




The integral is regularized by an UV cutoff A, and the UV divergences of the type 
A 3 log A and A0 2 , are absorbed into the total energy and mass, as usual. The poten- 
tial V e fr(0o) adds to the tree potential Vo = ^m^Q + iA^Q, and results in a first-order 
phase transition above T c at T\ = T c (l + £ 2 g 6 /187r 2 A), see Fig. |T[ 
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This conclusion has been questioned on various grounds. First of all, there is the 
experimental accessible smectic-A to nematic phase transition in liquid crystals. De 
Gennes' theory describing this transition can be mapped onto a Ginzburg-Landau 
model with the transverse fluctuations of the director field in the smectic-A phase 
playing the role of the gauge field ||. The transition is known to have a second-order 
regime 0. Furthermore, a lattice simulation carried out by Dasgupta and Halperin 
||, employing a duality map introduced by Peskin, Thomas and Stone f7|, showed 
that for small g-values the theory undergoes a second-order transition. 

To see the flaw in the argument leading to the previous conclusion, consider 
the Ginzburg region \tq\ < A 2 £q around T c in which amplitude fluctuations of the 
Higgs field can no longer be ignored. There, the London limit breaks down and 
the amplitude fluctuations of the Higgs field must be included. They have a mass 
term ^m^0 2 with mn = v^2|"^| the so-called Higgs mass. The ratio of the Ginzburg 
region to the interval T\ = T\/T c —1 in which the first-order jump takes place t g /ti ~ 
A 3 / q 6 ~ k 6 is proportional to the sixth power of the Ginzburg-Landau parameter k = 
m-ft/ \/2mA- As has been stressed by Kleinert ||, in type-II superconductors where 
k > 1/V2, the first-order jump takes place inside the Ginzburg region, where the 
amplitude fluctuations around <po are no longer smooth and of low energy, but driven 
by magnetic vortex loops. At the core of a vortex, the Higgs field vanishes making the 
effective potential (§) unreliable. As will be seen below, these fluctuations increase 
with increasing k. 

To describe a vortex loop, it is useful to first construct a vortex line with a 
source at some space-time point z, which acts like magnetic monopole. As such it is 
accompanied by a Dirac string. The electrodynamics in the presence of a monopole 
was first described by Dirac || who argued that F^ v = d^A v — d v A^ — F^ v plays the 
role of a physical field strength. The subtracted term 

/r dy\ 
dy x 5{x-y) = e^ A $ ds—5[x - y(s)], (3) 

with $ — 27r/g being the magnetic flux quantum, removes the field of the Dirac 
string running along some path L z from the location z of the monopole to infinity. 
The field F^ v (x) is a gauge field in its own right, with the corresponding gauge 



transformations changing the string configuration fll 



The monopole density is given by the divergence of the dual field F^ — ^€^ U \F^^, 
i.e., ^Fj(x) = §o5(x — z). In the Higgs phase, the Dirac string becomes the core 
of a vortex [|10j. 

The field equation obtained from the /1^-dependent terms in the Lagrange den- 
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sity, La = \F 2 + ^m\Aj l , is solved by 

d s yG(x-y)F^y) + %d,G(x-z); G(x) = J ^ + ^ , (4) 

with G(x) the massive propagator. Substituting the solution back into the Lagrange 
density, we obtain, ignoring the (diverging) monopole self-interaction, the action 
Sa = My f L ds of a relativistic point particle with mass 

My = imM / H^jj^r « >» - (5) 

In deriving this we employed the Lorentz invariance of the theory which allowed 
us to calculate the mass in the rest frame of the particle, and we took a vortex 
radius of the order of the coherence length . We see that in the Higgs phase, 

the vortex surrounding the Dirac string acts as the worldline of a point particle. 
Equation (f|) was first obtained by Abrikosov in the context of superconductivity 
|l"3|j . To see in which regime these excitations are important, we convert the mass 



My into a magnetic field via H Cl = My/$o, and compare it with the critical field 
H c = $ \m ( f > \mA/'2V2Tf which is obtained by equating the tree potential Vo to \H^. 
This field is a measure of the condensation energy which in turn sets the energy 
scale. It physically denotes the external field at which a type-I system can no longer 
resist the magnetic pressure and reverts to the Coulomb phase characterized by a 
perfect penetration of the field. The ratio of the two fields H Cl /H c = \nn/ 
shows that for increasing n, vortices become easier to excite. 

The field V(L Z ) describing the vortex is easily inferred by noting that in the 
functional-integral approach, a given field configuration is weighed by a factor e~ s . 
The Fj- dependent terms then defines the vortex field |14j]: V(L Z ) = exp[J d 3 x(^F^ + 
The Poisson bracket {F 0i (t, x), Aj(t, y)} = — <%5(x — y), stating that F ( 



0/ 



is canonical conjugate to implies that the vortex field and, thus, the point par- 
ticle carries one flux quantum as charge: {V, $} = &oV. We shall henceforth refer 
to the point particles as fluxons. The symmetry generated by the magnetic flux 
$ = / d 2 xFo is a global U(l) symmetry which is to be distinguished from the local 
U(l) gauge symmetry. As has been pointed out by Kovner, Rosenstein and collab- 



orators Ul5| , although the flux symmetry is of topological nature, the corresponding 
current being trivially conserved when expressed in terms of the gauge field, it 
is apparently not devoid of physical relevance. The central point of the dual formu- 
lation of the Abelian Higgs model is to describe the fluxons by a field theory where 
the flux symmetry becomes an ordinary Noether symmetry with a current whose 
conservation follows upon invoking the field equations. 
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Let us proceed to calculate the correlation function (V(L Z ) V*(L Z )) where V*(L Z ) 
describes an antimonopole located at z and its accompanying Dirac string. In the 
Lorentz gauge, we have 



(V(L g )V*(L a )) = J VA^exp [-§ J d 3 x (P* + m\A 



(6) 



where now d^F^x) = p(x), with monopole density p(x) = &o[5(x — z) — 8(x — z)]. 
Equation @ can be evaluated directly by integrating out the gauge field. For 
later convenience, however, we first linearize the Maxwell term via a Hubbard- 
Stratonovich transformation so that —\Ff L — » — \h 2 ^ l + ih IJL F ll . The physical relevance 
of the auxiliary field becomes apparent when considering its field equation, = iF^, 
showing that it represents the dual field. In this case, where the dual field is described 
directly by and not via an unphysical gauge potential, the flux symmetry is no 
longer a topological symmetry. The correlation function becomes 



(V(L Z )V*(L 2 )} = J Vh, exp j-/ d 3 x 
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[iri2 id,h v -d v h,y + -h^+ih^ 



(7) 

showing that the magnetic vortex couples with a charge g = $o^a to the massive 
vector field (the factor $o is hidden in the definition of the dual field F^, while 
the mass factor arises after rescaling the vector field so that the kinetic term has the 
usual coefficient 1/4). The remaining integral in (|7|) is Gaussian and can be easily 
carried out to yield: 

(V(L Z )V*(L S )} = exp - 6 e -^' L -l. (8) 

\ Q | Lizz I / 

In the Higgs phase, where the fluxon mass is nonzero, the second factor domi- 
nates and the correlation function tends to zero in the limit where the monopole- 
antimonopole pair is taken very far apart, so that the length of the vortex connecting 
the pair |L 22 | tends to infinity. This physically represents monopole confinement. In 
the Coulomb phase, where itla and therefore the fluxon mass is zero, we have on the 
other hand (V(L Z )V*(L 2 )) — > 1. On account of the cluster property of correlation 
functions, this implies that the vortex field develops a vacuum expectation value, 
thus signaling a proliferation of magnetic vortices. It should be noted that it is the 
Coulomb phase and not the Higgs phase where V(L Z ) develops an expectation value. 
It follows that the vortex field can be used to distinguish the two phases, i.e., it is 
an order parameter — or better, a disorder parameter |16| . 



It is instructive to consider the limit rriA — * directly in (0). Because of its 
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diverging coefficient, the first term must vanish, i.e., we have = (9^7, and 

(V(L Z )V*(L S )) = |P7exp{|d 3 a;[-i(^7) 2 + ^p]}- (9) 

This shows that in the Coulomb phase, V(L Z ) is represented as a normal field 
V(L Z ) = exp[i$ 7(z)] and not via a singular field Apparently, the Dirac string 
has become irrelevant here. The massless scalar field appearing in @ is the Gold- 
stone field of the broken global U(l) flux symmetry. The corresponding current 
jfi — df/y is no longer trivially conserved, but rather has become a Noether current 
whose conservation follows upon employing the field equation d 2r j = 0. Physically, 
7 represents the massless photon of the Coulomb phase. Since in two space and 
one time dimension the photon has only one transverse direction and thus only one 
degree of freedom, it can be represented by a scalar field. 

We are now in a position to write down the dual theory which physically describes 
a grand canonical ensemble of fluctuating closed magnetic vortices, of arbitrary 
length and shape. We will refer to such an ensemble as a loop gas. A neutral loop 
gas can be described by a massive scalar field theory with a repulsive two-particle 
interaction it^l 4 , accounting for the steric repulsion between vortices. In our case, 
the vortices interact in addition with the massive vector field h^, with coupling 



constant g = $o m A- This gives as dual theory [16, 15, I] 

c = ^r( 9 A - duKf + hi + 1(^ - *$oM^| 2 + ™JM 2 + <#l 4 , (io) 

where the -0-field is a disorder field whose Feynman diagrams are direct spacetime 
pictures of the fluctuating vortex loops in the original system (ID) |TB]| . It is minimally 
coupled to the magnetic field h^. The dual theory describes the system directly in 
terms of physical degrees of freedom; there is no local gauge symmetry. The old mass 
result (D appears in the dual theory as a one-loop contribution to the self-energy of 
the scalar field ip. The interaction strength of is determined by the dimensionless 
ratio g 2 /rriA, which is the inverse of that of the gauge field, g 2 /w,a ~ it^a/q 2 - This 
is a common feature of mutually dual theories. It is one of the reasons why the dual 
theory is more easily dealt with than the original Higgs model, where there is no IR 
fixed point in the weak-coupling regime, leaving only the strong-coupling regime for 
a second-order phase transition. 

Before carrying out the renormalization program we note that in the dual theory 
no cubic term is generated when the vector field is integrated out || . More precisely, 
repeating the steps leading to the effective potential V e g(4>), Eq. (@), we obtain now 

VeffW = / 7^ Hk 2 + m 2 A + 2g 2 \ip\ 2 ) = - J- (m 2 A + 2g 2 \^ 2 . (11) 

Jk<A [iT^y D7T 
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This has a Tayler expansion in \ip\ 2 and, in contrast to the previous case, no cubic 
term is generated in the Higgs phase characterized by a non-zero vector field mass. 
Since g = 0o m A, the effective potential V e fi(ip) vanishes as we approach the critical 
point by letting m A — > 0. Thus the massive vector field decouples from the ■0- field 
representing the vortices, and the Coulomb-Higgs phase transition reduces to a pure 
|^| 4 -theory, which is known to undergo a second-order transition. 

We now carry out the renormalization program, and write the bare theory fllPf ) 
as a sum of the renormalized Lagrange density and counter terms 

5C = Z f \(d„ - igh^Zg/Z^ipl 2 - |(9 M - igh^\ 2 + {Z^m\ fi - mj,)|^| 2 

+u(Z u - 1)H 4 + \{Z h - l){d^h v - dvhrf + \{Z h m% - m 2 A )h% (12) 

where hp = Z h 1/2 h ^, g = Z g l Z^Z^ 2 g , ip = Z^ 1/2 ip , and u = Z~ x Z\u^. [All 
quantities appearing in (|1(]) have been given an index to indicate that they refer to 
bare quantities.] An explicit one-loop calculation |l] , which uses the fixed-dimension 
approach of Parisi flL7|| , revealed that Z g = Z^, so that the minimal coupling is 
preserved at the one-loop order, and also that the mass term did not receive a 

1 /2 

correction, so that m A = Z h m^. Both these statements remain true to all orders 
in perturbation theory on account of the Ward identity known from QED which 
also operates in the massive vector case. As a result, both the mass and coupling 
constant of the dual theory renormalize in the same way, and they are linked via 
the equation g = (2ir/q)mA, with the electric charge q a free parameter which does 
not renormalize in the dual formulation. 

The double role played by wia as both a mass parameter and a coupling constant 
is essential to our renormalization group analysis We assume, as suggested by 
the first-loop calculation, that the phase transition is driven by the proliferation of 
vortices as signaled by a non-zero expectation value of the ^-field. On account of 
the relation g = $o m A, we can cast the /3-function P^(g 2 ,u) := dg 2 /<91n(m^)| , with 
u := u/m^ and the index indicating that go and u$ are kept fixed, of the properly 
scaled coupling constant squared, g 2 := g 2 /m^, in the form 



(3^{g 2 ,u) = g" 1 



-1 + 



dHg 2 ) 
dhx{m^) 
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^ d\n(m% ) d\n(Z h ) d\n{m A ) 
dln(m^) d\n(mA) dln(m^) 



Since rn\ ~ tt&q ~ r, the second term at the right-hand side is the inverse 
of the correlation-length exponent u, which indicates how fast the renormalized 
mass m^p tends to zero, ~ r v . The last term contains the 7h(g 2 , w)-function, 
dln(Zh) /d\n(mA)\o which yields the anomalous dimension r/h of the vector field 
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Figure 2: The /3-function of the dual theory (left figure), showing that the origin is 
an IR fixed point. For comparison, the expected form of the /5-function of the Higgs 
model (right figure) is included. 

when evaluated at the fixed point. That is, 



fy{§ 2 ,u) = g 1 



1 



2 dln(m A ) 
-1 + - +Jh{9 ,« 



<91n(m</,) 



(13) 



Without the double role played by m^, the coefficient of the g 2 term in the A/,(g 2 , «)- 
function would be —1, implying that the origin g = is an UV fixed point. Here, 
however, the coefficient is — 1 + \ jv which is positive if v < 1. In 3D, where v ~ |, 
this is the case, meaning that the origin g 2 = is IR stable and the massive vector 
field decouples, see Fig. ^} Hence, we are left with the pure |^| 4 -theory, just as we 
saw by expanding (pTTf) . 

This situation differs sharply from that in the original formulation, where the 
coupling q to the electromagnetic gauge field is assumed to possess a perturbatively 
inaccessible IR fixed point away from the origin, see Fig. ^. The Gaussian fixed 
point g* = u* = is a tricritical point of the Higgs model whose existence and 
location was first established by Kleinert f|. At this point, where the second-order 
transition changes into a first-order one, the vortices are completely free. 

We conclude that the Coulomb-Higgs phase transition constitutes a new uni- 
versality class. Whereas the disorder field behaves as in a pure |^| 4 -theory, the 
exponents of the magnetic field retain their mean-field values. 

Up to recently, this problem was of mere academic interest because the tem- 
perature resolution needed to probe the Ginzburg region could not be realized ex- 
perimentally. However, in various high-T c materials, this region has become fairly 
large \tq\ ~ 10~ 2 , and experimentally accessible. Experiments just carried out by 
Lin et al. [[HI on YBaCuO seem to confirm our predictions. They clearly show the 



mean-field behavior of the magnetic penetration depth. Earlier measurements of the 
specific heat had established the XF-model-like behavior of the heat capacity |19| . 
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